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Abstract 



The one-dimensional motion of N particles in the field of many incoherent waves is 
revisited numerically. When the wave complex amplitudes are independent, with a 
gaussian distribution, the quasilinear approximation is found to always overestimate 
transport and to become accurate in the limit of infinite resonance overlap. 
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The motion of a particle in the field of many waves is a fundamental process 
in classical physics [71181120] . usually discussed in the first chapters of plasma 
physics textbooks. For modeling purposes, this process can be described by 
the hamiltonian model 



where the particle has position x and momentum v (normalizing mass to 1 
and writing A m = qE m /k m ). The wave field has M ^> 1 waves, with a smooth 
dispersion relation associating a wavenumber k m , a pulsation u m and a phase 
velocity v m = u m /k m to each wave - usually determined by fixed properties of 
the environment, such as the geometry of the domain where waves propagate 
(then the k m are discrete). The actual spectrum of the wavefield, given by 
the complex amplitudes A m e lipm , is more easily tuned by the experimenter or 
affected by simple changes in the environment, e.g. excited antennas. 

1 Supported in part by a CNRS temporary position during this work. 

Preprint submitted to Communications in Nonlinear Science and Numerical Simulations30 May 2008 



transport 





m=l 



The dynamical system approach to this problem would be to prescribe a single 
choice for each wave complex amplitude. However it would be quite exceptional 
to control all waves (though this is e.g. the assumption underlying the standard 
map with A m = A , (p m = ip for all waves - see [5] for a discussion). 

Because the waves have different frequencies and velocities, it is generally 
unrealistic to assume their phases to be correlated. Their intensities are more 
easily observed, but both in nature and in the laboratory the accumulation 
of statistics on waves often involves only their average power spectra, not the 
detailed intensity data for each measurement run. Besides, if the waves are 
excited by a noisy source, as occurs in the self-consistent dynamics of particles 
and waves PIT2"] . one cannot expect a systematic reproduction of individual 
wave data among several experiments (this is the very meaning of a "noisy 
source" ) . 

We thus assume here that these complex amplitudes are random data, and in- 
vestigate the statistics of the particle motion in the resulting time-dependent 
random field. This dynamics is a "stochastic acceleration problem" for a "pas- 
sive particle" in weak plasma turbulence [6|19f21| . and its understanding is a 
prerequisite to a proper analysis of the case where the particle motion feeds 
back on the wave evolution [T2] . 

This hamiltonian generates equations of motion 



x — v, 

v = k m A m sin(fc m x uj m t fm) ■ 

m 



(0.2) 
(0.3) 



A simple case for the dispersion relation is 

k m = ho, u m = 2vr(m - M/2)/T, (0.4) 

for some ko, T. Then, in the limit M — > oo, the equations of motion yield 
for A m e liPm = A (with real A ) the well-known standard map [3]. The case 
where phases (p m are independent random variables uniformly distributed on 
the circle [0, 27r] while A m = A is given was investigated notably by Cary, 
Escande, Verga and Benisti [THHfT2] and occurs in the context of the random 
phase approximation. 

An important observation [51TT4"] on the motion of a particle in the field (10.31) 
is locality in velocity : the evolution of the particle when it has velocity v 
depends only weakly on the waves with a Doppler-shifted frequency u m — k m v 
much larger than their trapping oscillation frequency 2k m ^A m . In particular, 
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for a two-wave system the resonance overlap parameter 



2^ + 2^ 

Sl,2 = I 77 77-7 (0.5) 

becomes unity when there exists a velocity u for which uo 2 /k 2 — u = 2y/~A 2 ' 
and u — u)\jh\ = 1\fA[ (with k\ > 0, k 2 > 0, co>2 > ^i)- For many waves 
with overlap parameters s ^> 1, the relevant phase velocity range for waves 
influencing the particle is a "resonance box", with a width scaling as A 2 / 3 

To the extent that the phases and amplitudes of the waves are independent 
random variables, it is tempting to approximate the acceleration (10.31) by a 



white noise, with amplitude a = yE (k^A^) = k m ^J]EA 2 n , where the relevant 
mode m is the one with the current particle velocity (E denotes the math- 
ematical expectation, or ensemble average with respect to wave amplitudes 
and phases) H] 

This is the core of quasilinear theory [T0|17|ll8 ! ; 22j. Mathematically, one then 
interprets f l0.2p - fl0.3l) as a stochastic differential equation, after the necessary 
reformulation to make sense of divergent series ; this is correct in the case m G 
Z for dispersion relation (10.41) with gaussian independent complex amplitudes 
such that EA^ = k$ 2 a 2 . The particle velocity then has a brownian evolution, 
so that for < t < t' < T 

E (v t - vq) {v v - vo) = 2D Qh mm(t, t') (0.6) 

with the quasilinear diffusion coefficient Dql = a 2 T/4. In particular, 

EA^ = yt (0.7) 

for Avt = Vt — vq. However, the particle evolution for t > T may show a strong 
correlation to its motion for < t < T because the waves are periodic in 
time [TT1fT3] . and the dispersion relation (10 .4p may generate a strong spatial 
correlation - although it just accounts for the fact that in the limit of a dense 
spectrum (with a smooth dispersion relation) the relevant waves acting on a 
particle at any time also have almost the same wavelength^] 



2 Phases do not appear in a (nor in s) because the relative phase of two waves 
ip m + u> m t — tp n — uj n t varies uniformly over time (hence ip m — (p n can be absorbed 
in the choice of the time origin) . 

3 There is a large body of literature on the case of incoherent waves with no dis- 
persion relation. Then the sum becomes a double sum n and one varies 
wavenumbers k n independently from pulsations to m . This space-time stochastic en- 
vironment is more noisy that our model and may also be considered to motivate a 
quasilinear approximation. 
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In this note we compare numerical properties of the velocity evolution for fixed 
amplitude with the case of random wave amplitudes : let A m e 1(fim = C m + iS m . 
In the latter case, the amplitudes are drawn independently in such a way that 
their squares have exponential distributions, with equal expectation Ei^ = 
k§ 2 o 2 . Along with uniform distribution of the phases, this is equivalent to 
gaussian distribution of the complex amplitudes. Indeed, the probability den- 
sity fp for (p is uniform, and the exponential law for A 2 yields for A the density 
at a > 

(a < A < a + da) _ P(a 2 < A 2 < a 2 + 2a da) _ ^_ a 2 /Ej4 2 2a 



Ua) ~ d^ " d^ ~ e ~ ea*- 

(0.8) 

The probability density for (C, S) is thus (noting that a 2 = c 2 + s 2 and dcds = 
a da dip) 

_ f A (a)f ip ( ip )dad V _ 2ae~ a2 / EA2 dad<p _ e -(^ 2 ^ 2 )/^ 2 
/cMC ' Sj dcds 2ttEA 2 dcds vrEA 2 1 j 

which is the density for a gaussian vector (C, S) with zero expectation and 
■EA 2 /2 



covariance matrix 



EA 2 /2 



For the simulations, we let ko = 1 and T = 2tt so that o = s 2 /16. Given the 
value of a, or equivalently of s, we draw at random iV p h a = 400 sets of wave 
data. For each set of wave data, we follow the evolution of iV or b = 50 particles, 
released initially at random points (x,p) in the strip < x < 2tt, —0.5 < 
p < 0.5. Trajectories are computed with a reversible symplectic integrator 
up to a time, large enough for the motion to exhibit possible departure from 
quasilinear statistics, but short enough to ensure that the particles remain far 
from the boundaries of the wave spectrum, which have velocity ±M/2 (where 
we take M up to 800) : for s = 13, the boundaries are more than 3 standard 
deviations (which is about 120 as observed in Fig. [1]) away from the initial 
particle velocity. Some trajectories are computed backward from their final 
point to control numerical accuracy : the accumulated error is about 10~ 5 
for s = 3.5 for the runs of Fig. CD and El though it deteriorates rapidly with 
increasing s. Our calculations for the equal A case reproduce the findings of 
Refs [TB] . 

Statistical averages discussed below (denoted by (■)) are performed over all 
particle and wave data for the same s value. We focus on the second and fourth 
moments of particle velocity as functions of time. To the numerical average 
92(t) = (Af 2 ) we fit a linear approximation D + 2D c fft, which defines an 
effective diffusion coefficient. Such a plot is displayed in Fig. [H The quasilinear 
prediction 2Dq L t is a straight line. The numerical average g 2 is plotted, along 
with two lines estimating one standard deviation on either side, viz. g 2 ±(Av 4 — 
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A = 1 0.56, M = 800, s = 1 3.00, dt = 0.000063, npha = 400, norb = 50 
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Fig. 1. Dependence of (Av 2 ) on time for random amplitudes (Gauss) and for equal 
amplitudes (vl=cst). 

<?f) 1//2 . The narrowness of the channel so constructed indicates the statistical 
accuracy of the plotted function g^. 

Specifically, Fig. [T] displays our results for s = 3.5 and for s = 13. The non- 
overlap regime (s < 1) allows no large scale transport in velocity, because 
the particle phase space contains invariant Kolmogorov-Arnol'd-Moser tori. 
Indeed the case s = is integrable (it corresponds to free particle motion). 
For increasing overlap parameter value, one expects the transport to become 
increasingly chaotic, and indeed the particle dynamics typically has a Liapunov 
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Fig. 2. Ratio -D e ff/-f ql as function of overlap parameter s for random amplitudes 
(Gauss) and for equal amplitudes (A=cst). The lines are guides to the eye. 

exponent scaling like A 1 / 3 ~ s 2 / 3 for large s [12J. One may then expect the 
effective diffusion coefficient to increase towards the quasilinear estimate, -Dql, 
which corresponds to a pure stochastic particle behaviour. However it has 
been found [TflfTiZ] that for the equal amplitude, independent phases wave 
data, the particle velocity undergoes an enhanced transport, with an effective 
diffusion coefficient up to 2.5 times the quasilinear estimate, and the latter 
becomes accurate only in the large s limit. It must be noted that this enhanced 
transport shows up only after time T = 2tt, i.e. only on a time scale for 
which the time periodicity due to the discrete nature of the wave spectrum 
significantly affects the particle motion - actually the motion for short times 
is essentially quasilinear thanks to the wave stochasticity [T2TT6] . We also 
checked that the force correlation function on the particle is essentially zero, 
except for times near integer multiples of the period T. 

Our simulations show that the enhanced transport is not observed for inde- 
pendent wave amplitudes. In the latter case, the velocity variance grows with 
time at a rate never exceeding the quasilinear value. As shown in Fig. [2] this 
observation holds for the whole range of overlap parameter values sampled, 
in agreement with the naive prediction that "increasing chaos" should make 
dynamics more similar to "pure noise" . 

The noisiness of the particle evolution is also assessed using a higher-moment 
test. If the velocity distribution at time t were gaussian, the ratio (At> 4 )/(3g|) 
would equal 1. Fig. [3] displays numerical evidence that for the equal wave 
amplitude case this ratio remains rather below 1 (the velocity distribution has 
smaller tails than the gaussian with equal variance), whereas for the gaussian 
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Fig. 3. Ratio (Av 4 )/(3(Av 2 ) 2 ) versus time for random amplitudes (Gauss) and for 
equal amplitudes (A=cst). 

wave data it is somewhat above 1 (the distribution has stronger tails than the 
gaussian with same variance). 

For both wave data types the limiting behaviour for s — > oo is quasilinear, in 
agreement with theoretical [T|12|15j and mathematical [Tl~1fT3] arguments. 

These results confirm that random phases only are not sufficient to substanti- 
ate the quasilinear approximation for the stochastic acceleration problem. A 
gaussian wave spectrum may seem closer to the ideal view of white noise for 
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one period T, but over longer times correlations also build up, driving trans- 
port away from the quasilinear approximation, especially when the overlap 
parameter has moderate typical values (1 < s < 5). 

Many fruitful discussions with F. Doveil and D.F. Escande are acknowledged. 
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